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A validation of the matrix dissipation option in OVERFLOW was completed, and a 
single set of constants was found which yield good inviscid shock capturing, viscous heat 
and skin friction predictions, and work over a wide Mach and Reynolds number range. 

This scheme works well with multigrid, and produces fast convergence to steady state for 
2 and 3 dimensional problems, with accuracy comparable to Roe upwinding, and without 
’carbuncle’ problems for the blunt body flowfields. 


Introduction 

OVERFLOW 1 ’ 2 has been mainly applied to flow 
past subsonic, transonic and supersonic complex ge- 
ometry vehicles. An investigation of the flow solver in 
the ideal gas regime is the first step in extending the 
capability of OVERFLOW to include non-equilibrium 
chemistry flowfields. The requirements for the flow 
solver were 

• capable of handling strong shocks, yielding good 
post-shock conditions 

• capable of providing a starting solution 

• good skin friction and heat transfer predictions 

• extensible to non-equilibrium chemistry flowfields 

Matrix dissipation 3 mimics a TVD upwind biased 
scheme, within an underlying central spatial differenc- 
ing framework. A commonly cited problem is the trial 
and error procedure needed to choose values for the 
smoothing parameters to ensure stability, and obtain 
solutions of high accuracy. A set of universal pa- 
rameters has been found which will provide sufficient 
stability to converge inviscid blunt body flowfields, 
yielding correct shock strength for even moderately 
coarse grids. These parameters eliminate ‘carbuncle 1 
problems in these flowfields, and still provide accu- 
rate predictions of skin friction and heat transfer for 
flat plates for subsonic to supersonic flowfields. A 
further advantage is that these parameters, and the 
central /matrix dissipation scheme work well with the 
existing multigrid framework within OVERFLOW to 
rapidly converge to steady state. The scheme has been 
tested on moderate to high hypersonic Mach number 
flowfields, on 2D and 3D problems. 
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Method 

The scheme chosen here combines 2nd order central 
pacial differencing with the Swanson /Turkel matrix 
dissipation. When using this method, 4 adjustable pa- 
rameters must be chosen: the coefficients of the 2nd 
and 4th smoothing^ and K4 in ref, 3 DIS2, and DIS4 
in the OVERFLOW nomenclature), and two eigen- 
value limiters(V €ri , V €a in, 3 VEPSL and VEPSN for 
OVERFLOW). A single set of constants was found 
which satisfied the requirements enumerated above. 
Roe averaging was used to obtain the matrices used 
in the dissipation scheme'since we are attempting to 
obtain a shock capturing scheme for high normal Mach 
numbers. 

The constants used to obtain the results in this pa- 
per are: 

DIS2 = 2 
DIS4 = 0.1 
V Ca = 0.3 
Vei =0.3 

The P ulliam- Chausee diagonal scheme, coupled 
with multigrid and grid sequencing is used for a relax- 
ation scheme, providing fast and robust convergence. 

Starting flowfields with detached shocks is the only 
time when the parameters of this scheme require ad- 
justment. To obtain a starting solution, the flowfield is 
started from a uniform freestream on a very coarse grid 
(15 of every 16 points in the fine grid are dropped) with 
both eigenvalue limiters (V €l , V €ti ) set to unity. Set- 
ting the limiters to unity makes the dissipation mimic 
a scalar dissipation scheme. This produces a flowfield 
with a smeared, but correctly located shock, and does 
so very cheaply. The residuals should be dropped at 
least 3 orders of magnitude in this step. The solution 
thus obtained is quite amenable to allowing a shock 
fitted grid to be obtained with grid adaptation, if de- 
sired. The results of this paper were obtained without 
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a) Pressure b) Temperature c) Density 
Fig. 1 Inviacid Moo - 32 Blunt Body Flowfield 


a scalar dissipation scheme* This produces a flowfield 
with a smeared, but correctly located shock, and does 
so very cheaply. The residuals should be dropped at 
least 3 orders of magnitude in this step. The solution 
thus obtained is quite amenable to allowing a shock 
fitted grid to be obtained with grid adaptation, if de- 
sired. The results of this paper were obtained without 
any grid adaptation. After the starting solution is 
obtained, the solution is restarted with both V € pa- 
rameters set to the recommended value of 0,3. Grid 
sequencing is recommended as well as multigrid. 

Results 

Inviacid Circular Cylinder 

The overall flowfield is shown in Fig.l, with contours 
of the log of pressure. 

This case is the inviscid flow past a blunt cylinder 
at various Mach numbers. This case was run at 3 
Mach numbers: 8, 16 and 32, to obtain solutions with 
moderate to extreme shock pressure jumps. A 101 x 
101 grid was the finest grid run, and cases with 51 x51 
and 26 x 26 are also shown, to exhibit the schemes 
ability to handle coarse grids robustly. 

Convergence 

The M = 32 case is considered as an example here, 
as it is the most challenging case to converge. An 
initial run to set up the starting solution was done as 
described above, with 600 iterations serving to lower 
the residual by 4 orders of magnitude. 

Next the solution was restarted with the limiter pa- 
rameters reset to the recommended values of V <n = 
V Cl = 03, and through 3 stages of grid cycling which 
keep the “depth” of the multigrid method at 2. All 



these cycles were run using constant CFL of 0.8, and 
resulted in a fine grid residual reduction of 5 orders 
of magnitude for the first 1000 fine grid iterations. 
The residual histories are plotted in Flg.2 with resid- 
uals shown as a function of cpu seconds. The relative 
cheapness of the coarse grid cycles is graphically evi- 
dent. 

Fig.2 shows the ability of multigrid to accelerate 
convergence for this case. The residual history as a 
function of CPU seconds is shown for the fine grid so- 
lution with and without multigrid. Multigrid seems 
to provide faster and more robust convergence for the 
other flowfields we have investigated as well. That the 
multigrid scheme actually worked more reliably and at 
higher CFL numbers than non-multigrid for an invis- 
cid flowfield with strong shocks was unexpected, but 
welcome. For viscous flowfields, the multigrid scheme 
greatly accelerates convergence. 

Smoothing Parameters 

The choice of smoothing parameters was fixed 
largely by this problem, with the “carbuncle prob- 
lem” as the determining factor. When the 4th order 
smoothing value (DIS4) is lowered by a factor of two 
(from .1 to .05;Fig.3), or with the linear ” limi ter” set 
to V Cl < 03, the temperature, density velocity and 
contours of this case begin to show signs of the car- 
buncle problem, where insufficient smoothing allows 
the solution to develop slope discontinuities in deriva- 
tives post shock in the stagnation region. 

The 2nd order smoothing parameter, DIS2, can be 
increased from its recommended value. The only ef- 
fects of increasing DIS2 are that the shock is smeared 
over more points, and the stagnation enthalpy (which 
should be constant for this problem) is less well con- 
served, as seen in Fig 4). Fbr example, stagnation 
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the solution to develop slope discontinutities in deriva- 
tives post shock in the stagnation region. 

The 2nd order smoothing parameter, DIS2, can be 
increased from its recommended value. The only ef- 
fects of increasing DIS2 are that the shock is smeared 
over more points, and the stagnation enthalpy (which 
should be constant for this problem) is less well con- 
served, as seen in Fig 4). For example, stagnation 
enthalpy and stagnation pressures are altered by less 
than 1% if DIS2 is increased by a factor of 5. If DIS2 is 
decreased by a factor of 2, the scheme fails to converge. 



A similar dependence on the nonlinear eigenvalue 
limiter V €n < 0.3 parameter is also shown. It is pos- 
sible to decrease V £n to 0.2, without creating slope 
discontinuities in the stagnation region, but conver- 
gence is not as rapid, with no obvious improvement 
in solution accuracy (Fig 4)). Increasing any of the 
smoothing parameters above their recommended val- 
ues actually slows the convergence process. Although 
these results have been shown for the M = 8 case, 
these observations also hold true for the higher Mach 
numbers, at least up to M = 32. 

Accuracy and Grid Requirements 

In these discussions, the static pressure and stagna- 
tion enthalpy will be focused on, as these two variables 
will determine the wall boundary layer evolution for 
boundary layer flows. The accuracy to which these 
variables can be predicted directly affects the accuracy 
of the skin friction and heating predictions for viscous 
flows. 

Stagnation line profiles of pressure and stagnation 
enthalpy for M = 32 case for solutions on coarse, 
medium and fine grids are shown in Fig 5. This case 



has extremely strong shocks, with a shock pressure 
rise of 3 orders of magnitude. The stagnation point 
pressure is predicted to within 1% on all grids, with 
the finest grid within 0.1% of the theoretical value. 
Even on the coarse grids, the scheme produces good 
post-shock and stagnation predictions. The error in 
stagnation enthalpy is less than 2% on the finest grid 
at the stagnation point, and for the fine grid, the error 
is less than 0.2%. This ability to produce reasonable 
solutions on coarse grids is presumably the reason for 
matrix dissipation’s good performance when coupled 
with multigrid. 

It should be emphasized that all these solutions were 
obtained with the same values of DIS2, DIS4, V €n 
and V Cl . The question now is whether with these 
smoothing parameters, viscous boundary layers can be 
accurately predicted. Thus, we are led to the flat plate 
test cases, where the ablity to accurately predict the 
boundary layer flow may be assessed. 

Compressible Laminar Flat Plate 

For this case, without a detached shock present in 
the ftowfield, it is no longer necessary to obtain a start- 
ing solution (cycle 1 in the sequence outlined above). 
A single run of OVERFLOW is sufficient to obtain 
the solution for this case, and is made with the rec- 
ommended parameters, with the full multigrid (grid 
sequencing) of 200 iterations on a coarse grid, 400 on 
a medium grid, and 400 iterations on the finest grid. 
Convergence was checked by running 1000 iterations 
further on the fine grid. 

The laminar flat plate was run at M = 2, and 
M = 4, with adiabatic and constant temperature wall 
conditions, and a freestream temperature of 51 7.8° R. 
The grid for this case was 129 by 129. The skin friction 
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A single run of OVERFLOW is sufficient to obtain 
the solution for this case, and is made with the rec- 
ommended parameters, with the full multigrid (grid 
sequencing) of 200 iterations on a coarse grid, 400 on 
a medium grid, and 400 iterations on the finest grid. 
Convergence was checked by running 1000 iterations 
further on the fine grid. 

The la m i n ar flat plate was run at M = 2, and 
M = 4, with adiabatic and constant temperature wall 
conditions, and a freestream temperature of 51 7.8° R. 
The grid for this case was 129 by 129. The skin friction 
comparisons for the T w = Too cases are shown In Fig. 6 
and Fig.7. The wall heating rate (Ch) comparisons 
are shown in Fig.8 and Fig.9. The heating rate pre- 
dictions and skin friction are In very good agreement 



with theory over the entire plate. These figures also 
demonstrate the independence of these predictions to 
variations in the 2nd order smoothing coefficient, as 
the plots show open symbols with the recommended 
smoothing values, and filled symbols solution results 
with DIS2 increased by a factor of 2. 

Now, the inviscid shock capturing ability and the 
ability to model viscous flows has been demonstrated 
in isolated cases. We are now led to a simple case 
involving both shocks and viscous effects, the viscous 
circular cylinder. 

Viscous Circular Cylinder 

Given this good agreement for the inviscid blunt 
body and viscous flat plate, a natural question is how 
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Fig. 10 Stagnation Line Profiles for Viscous Circu- 
lar Cylinder, Moo = 11.436, Tw — 1200. K,Too — 265 K 


51 x 51 grid was obtained from the original grid by 
removing every other point from the finer grid, and the 
26 x 26 grid was obtained similarly from the 51 x 51 
grid. 

The stagnation line profiles for this case are shown 
in Fig. 10. The adiabatic wall temperature for this case 
is 6146.53°K, so that if the freestream enthalpy is 1% 
low downstream of the shock, the heat transfer predic- 
tions would be expected to be 1.3% low. The fine and 
medium grids have an enthalpy about 2% low, and 
the coarse grid has an total enthalpy 10% low. All 
grids predict the same pressure profile, but the post 
shock temperature is lower and the density higher on 
the coarse grid, consistent with the low total enthalpy 
post shock. 

The stagnation point heating rates are predicted 
reasonably well on all grids (Fig. 11), with the fine grid 
l.% low, the medium grid 1.1% high and the coarse 
grid 10% high. The overprediction of heat transfer on 
the medium and coarse grids is due to the coarser wall 
resolutions and the slight overshoots in total enthalpy 
at the boundary layer edge. 

The fine grid in this case required a higher 2nd or- 
der smoothing (DIS2=4, twice the “recommended” 
value) in order to achieve convergence in c h - The 
heating rate prediction was far more sensitive to so- 
lution convergence than quantities like pressure were. 
Small variations in the solution which could not be 
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Fig. 11 Heating Rate Predictions, for Viscous 
Circular Cylinder, Moo = 11.436, T w = 1200. K,Too = 
265° K 




Fig. 12 Cylinder Flare Flowfield, Moo — 9.91,T W — 
293°K, Too = 51°K 

discerned to plotting accuracy in the stagnation line 
profiles could make large variations in the heat trans- 
fer coefficient. 

Now, the ability to handle attached viscous flows 
with shocks has been demonstrated. The next element 
of flow complexity to be checked is the ability to model 
high speed separated flows. 

Viscous Cylinder-Flare 

This case 4,5 is a hollow, finite length cylinder, with 
a 10° flare. The flowfield exhibits a laminar separa- 
tion bubble. This flowfield requires a large number 
of iterations to converge as the laminar separation 
slowly creeps upstream. Fig. 12 shows the geometry 
along with the pressure field contours for a converged 
321 x 301 grid solution. 

A grid convergence study was completed for this 
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Fig. 12 Cylinder Flare Flowfleld, Moo * 9 . 9 1 , T w » 
2?3°K,Too-5TK, p«> -52Pa 



most clearly shown in Fig.14, which also shows a simi- 
lar grid convergence. The heat transfer predictions are 
similarly grid converged(Fig.l5). 

As noted previously, this flowfield required many 
iterations to converge, with the laminar separation 
moving slowly up the cylinder in the subsonic portion 
of the boundary layer. Over 10,000 iterations for the 
321 x 301 grid, and 6,000 iterations for the 161 x 151 
grid were required to reach the steady state solution. 
The solutions were run out (to 40,000 iterations for 
the 321 x 301 solution) to confirm the convergence of 
the relaxation process. Sage 7 was used to obtain the 
621 x601 grid, which also gave a very respectable start- 
ing solution. This solution was run out an additional 
10,000 iterations to confirm its relaxation convergence. 

The pressures predicted by both codes are above 
the experimental data upstream of the flare, simi- 
lar to the experimental/computational comparisons by 
Chanetz. 5 The pressures predicted on the flare are 
in good agreement with the data, as well as each 
other. The separation location predicted by both 
codes (X/L=0.7) is slightly upstream of the reported 
experimental value(X/L = 0.76 ± 0.01), but the reat- 
tachment point (X/L=1.35) agrees with the experi- 
mental value (X/L = 1.34 ± 0.015), The heat transfer 



Fig. 14 Cylinder Flare Skin Friction 



predictions are In good agreement with each other, and 
experiment. 

With the ability to handle this level of flow com- 
plexity in two dimensions, the next logical step is 3 
dimensional flowfields. 

3D Flowfields' 

The scheme described here has been successfully ap- 
plied to X38 forebody, “Sharp”, Shuttle and X33 grid 
systems. Here we will focus on the solution about 
the X33 geometry, as this used the complex geome- 
try (chimera) capabilities of OVERFLOW. 

The starting point for this grid system was the 
patched grid used to solve the flowfield with GASP. 
Preliminary computations utilized a chimera grid sys- 
tem obtained by simply extending the patched system 
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patched grid used to solve the flowfield with GASP. 
Preliminary computations utilized a chimera grid sys- 
tem obtained by simply extending the patched system 
grids into each other to provide sufficient overlap for 
the chimera scheme to work. 

Solutions obtained on this grid system appeared ad- 
equate, with pressure and heat transfer plausible over 
the vehicle surface. There was, however an under- 
lying problem, in that the chimera overlap surface 
included the stagnation point flow. With the minimal 
overlap created with this grid system, the surface pres- 
sure predictions were very good, but the heat transfer 
predictions, even when the grid appeared to be fully 
converged, were off by a factor of two. 

A new nose grid fully enclosing the stagnation point 
was rapidly generated using the chimera grid tools 
package (Fig. 16) and the solution was re-run. The 
pressure predictions obtained agreed with those ob- 
tained on the first grid system, and the heat transfer 
predictions came into line with those expected from 
previous GASP solutions, 7 8 and experiment. 




This solution required a residual reduction of 8 
orders of magnitude in order to provide good heat 
transfer estimates. Multigrid was again quite useful, 
although it had to be shut off in the aft grid which con- 
tains the wings and winglets. The upstream chimera 
join (near the nose) is a true chimera interface, where 
there are no coincident points. The downstream inter- 
face was created simply to break the body grid into two 
zones, in order to facilitate the parallelism of the solu- 
tion process. This interface has coincident surfaces. 

The surface flow patterns are continuous across the 
chimera zone boundaries along the surface, as can be 
seen in Fig. 19. Furthermore, the zonal boundaries in- 
tersect the bow shock on all three grids (Fig.20). On 
both of these boundaries, the chimera system provides 
sufficient solution continuity even for heat transfer es- 
timates. 


Fig. 18 X33 Centerline Heat Transfer, Moo - 5.98, 

a = 40°, Re L =3.3 x 10 6 

Discussion 

This work extends a number of innovative exten- 
sions of subsonic technology into the high speed world. 
Matrix dissipation has been used successfully in tran- 
sonic flowfields, and can yield predictions which are 
of the same quality as upwind MUSCL schemes, at 
much lower computational cost. This was found to 
be true even at high Mach numbers, with M 32 
normal shocks nicely modelled, even on coarse grid 
systems. Furthermore, one set of constants was found 
which can reliably used on all the grid systems tested 
here. While some attention was paid to keeping grid 
stretching below 1.2, and reasonably orthogonal grids, 
the grids used were generated with standard grid gen- 
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Fig. 19 X33 Surface Pressure Contours, Moo — 

5.98, a *40% Re L =3Jx 10* 



Fig. 20 X33 Symmetry Plane Pressure Contours, 

Moo - 5.98, a * 40% Rc L - 3 3 x 10* 


high speed flows to speed time to solution. Multigrid 
as part of the solution process does not seem to be 
as widespread, especially in the structured grid world. 
Multigrid was found to be valuable in speeding up the 
convergence process, on both Euler and viscous grid 
systems. Fbr solutions where extremely large residual 
reductions are required, such as heating rate predic- 
tion cases, this is invaluable. 

Another piece of transonic technology which worked 
well at hypersonic speeds was the Pulllam-Chausee 
diagonal algorithm. Although coupling this solution 
algorithm with chemistry is problematic, for solu- 
tions where the mean flow and chemistry systems can 
be solved in an uncoupled manner, this algorithm 
promises large computational cost savings. 

The solution of the X33 system shows the feasibility 
of using the chimera system for high speed flowfields. 
The ability to use chimera to model complex geome- 
tries greatly facilitates the ability to obtain predictions 
for flight vehicles, and allows the tailoring of grid sys- 
tems component by component. For rapid evaluation 


of an evolving system design, this capability is ex- 
tremely important. 

Conclusions 

Matrix dissipation combined with multigrid and grid 
sequencing shows promise as a method to solve hyper- 
sonic Mach number flowfields. The solutions obtained 
are as accurate as upwinding schemes, and appear 
to converge more quickly, especially when combined 
with multigrid. The chimera methodology also shows 
promise in this regard. In addition to providing a sim- 
ple means to model flows about complex geometries, 
the solutions obtained can be made sufficiently contin- 
uous across zonal boundaries. 

This work is part of a overall effort to use OVER- 
FLOW to predict vehicle re-entry flowfields. The 
scheme provides a fast, robust, and accurate means 
of solving the Navier-Stokes equations at hypersonic 
Mach numbers. 
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